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Elastic scattering in a quantum wire has several novel features not seen in ID, 2D or 3D. In this 
work we consider a single channel quantum wire as its application is inevitable in making devices 
based on quantum interference effects. We consider a point defect or a single delta function impurity 
in such a wire and show how some of these novel features affect Friedel-sum-rule (FSR) in a way, 
that is quite unlike in ID, 2D and 3D. 
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The density of states in a mesoscopic sample is very 
important to understand mesoscopic transport phenom- 
ena and also thermodynamic properties. It is believed 
that the properties of a mesoscopic sample, connected to 
leads, can be formalized in terms of the scattering ma- 
trix Ola- The Friedel -sum-rule (FSR) relates the den- 
sity of states (DOS) to the scattering matrix and can be 
stated as 6 f {E 2 ) - f (Ex) w nN(E 2 ,E 1 ) Q. In ID, 2D 
and 3D, the equality being approximate is almost exact 
in the WKB regime where generally transport occurs. 
Here N{E 2 ,E\) is the variation in the number of states 
in the energy interval [E% , E2] due to the scatterer and 
f (E) = IX^-fj = jiln(det[S]) 0. S is the scattering 
matrix and e l ^ j , j=l,2,....n are the n eigenvalues of the 
unitary matrix S. In differential form the FSR can also 
be stated as 

^f( E ) = ^ l ™( det iS})*K[p(E)-p (E)} , (1) 

where [p(E) — po(E)] is the variation of the DOS or the 
difference in the DOS due to the presence of the scatterer. 
p(E) and po(E) can be found by integrating the local 
DOS (LDOS) p\x, y, z, E) and p (x, y, z, E). 

In 1994, Ref.yJ explored a relation between LDOS and 
the scattering matrix. We refer to it as the Biittiker- 
Thomas-Pretre (BTP) sum rule which could be exact 
but is of limited practical value as far as the global DOS 
is concerned, because of the following reasons. First of 
all if one wants to integrate the LDOS to find the DOS, 
then it is very cumbersome and may not be at all pos- 
sible for complicated potentials. Secondly, in this local 
formalism, one has to take the derivatives of the S matrix 
elements with respect to the local potential and in some 
cases one encounters problems in including the non-local 
effects and non-local disturbances in the LDOS. Thirdly, 
if one wants to find the change in the S matrix due to 
an infinitesimal change in the local potential, then one 
has to know the change in the local wave-function due 
to that infinitesimal change in the local potential and 
then use the transfer matrix multiplication method or 
any other equivalent method to find the change in the S 
matrix. So then integrating the LDOS to find the DOS 
is just equivalent to integrating the wavefunction to find 
the DOS. Another such relation between the LDOS and 
the scattering matrix is given by Titov and Schomerus 



Whereas, the power of the FSR lies in the fact that 
it gives this integrated quantity straightaway, without 
having to tamper the internal structure of the potential, 
but just by infinitesimally changing the energy of an ex- 
ternal probing particle. Although, the BTP sum rule is 
of use in certain cases, it is definitely not a competitor of 
or a substitute for the FSR. No doubt, even the authors 
of Ref. [1], much after the proposition of the BTP sum 
rule in 1994, refers to Eq. 1 as the FSR, and are worried 
about to what extent it is valid 0,0- 

The scattering matrix of an impurity in a quantum 
wire have very unusual features that were not realized 
until very recently 0, . In a single channel quantum 
wire, in the presence of a single attractive impurity, taken 
as a negative delta-functionpotential, the transmission 
probability can go to zero [8j for some finite energy of 
the incident electron. At the corresponding energy, the 
scattering phase-shift shows a discontinuous jump(slip) 
by I" 13- It was shown that in the single channel case 
the Friedel phase Of is not affected by the discontinuous 
phase drops [rj[To|. 

In the multichannel case, when the unitarity of a par- 
ticular channel is not present and the electron can es- 
cape to a different channel, the transmission zeroes are 
replaced by minima and the discontinuous phase slip 
by 7T are replaced by continuous and less than 7r phase 
drops ^lj- K was also shown that in the multi channel 
case too the Friedel phase Of is not affected by these con- 
tinuous phase drops [Tlj. However, dS L^ may not bear 
any resemblance to n[p(E) — po(E)] . 

In this work we intend to study the FSR (Eq. 0J) for a 
single channel quantum wire with a delta-function impu- 
rity. The single channel case being the most important 
because it is in this regime that one can really control 
the quantum interference effects and use them to build 
mesoscopic devices|l2j and point impurities like the delta 
function potential are always present. Besides, to study 
the transport across a quantum dot connected to two 
ideal leads on two sides, most theoretical works model 
the dot by a single bound state at the site of the dot 
as the Coulomb blockade makes the other levels of the 
dot to be very far away. As an attractive delta poten- 
tial is capable of creating such a single bound state, it 
was used in Ref. to explain the Fano resonances in 
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quantum dots and the unusual features of the scattering 
phase shift observed across the quantum dot. Besides, 
the present study provides general understandings, and 
we also understand the system and the results of Ref. 
better. 

Scattering due to a spherical defect leads to asymp- 
totic wave-functions that are plane waves and whose ra- 
dial part is of the form V'fc,;( r ) ~ ~ Sin(kr — + rji) 
|13| . Here, k is the wave- number, I is the angular mo- 
mentum quantum number and rji is the phase shift due to 
scattering. According to Friedel [I3> to approximately 
count the number of states created by the impurity, con- 
sider a large sphere of radius R with a defect at the cen- 
ter, and impose a condition ipk,i(R) = 0. Thus we obtain 
kR — y + Vi = n7T ■ The states k thus obtained are real 
and not complex as it should be. Generally, scattering 
states have open boundary conditions that lead to their 
characterization in terms of complex energies. The imag- 
inary part of their complex energy is called self energy 
and it arises because of the fact that the states can leak 
out to infinity and get absorbed by some detector there 
irrespective of boundary conditions ^3 • I n this case also 
the interaction of the states within the large sphere of 
radius R with the region outside the sphere will lead to 
a self energy. In absence of the scatterer, the scattering 
phase shift rji = and thus one gets Eq. [f| . 

A more rigorous derivation, including the self energy is 
given by Buttiker et al • Earlier treatment of Dashen, 
Ma and Bernstein 01 an d that of Avishai and Band [l^ 
are also reviewe d by these authors, in some of their works 
like Refs. 0, 0, Il6j . These treatments are only valid for 
large systems. Souma and Suzuki [13 provides a straight 
forward extension of the work by Avishai and Band to 
quantum wires and also suffers from the same drawback. 
Ref. gives 

g^e f (E) + Im Tr G a — = Im Tr G a 



n[p(E)-p (E)} . 



Here G a 



E-H 



system 



t a (E) 



(2) 



is the ad- 



vanced Green's function and E° is the corresponding self 
energy. Note that G° is the advanced Greens function 
for the system alone, where the modifications in the sys- 
tem due to the presence of the leads is included. Hence 
apart from this Im TrG a which is equal to the inte- 
grated disturbance in the LDOS created by the impurity, 
i.e., ir[p(E) — po(E)], there will be some disturbance in 
the LDOS in the leads which will depend on the Greens 
function of the lead and how it is affected by the sys- 
tem. This contribution is not important and also it gets 
screened away very easily as it is very small and the leads 
being ideal, carrier concentration is very high in the leads 
(this is often referred to as non-polarisable leads). The 
only assumption required to get Eq. JTJ is to neglect 



i.e. the energy dependence of the self energy. Now 0] 



dE^ ~ Jfr. [ t p9p t p] 
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t p is the matrix element which couples the leads to the 
sample. In ID, 2D and 3D r p is independent of energy 
in the WKB regime. g p a being the local Greens function 
of a semi-infinite ideal wire, its energy dependence is not 
very important and thus in ID, 2D and 3D, in the WKB 
regime, FSR [QJ] is almost exact. There will be viola- 
tions in the non-WKB regime, but since transport occurs 
in WKB regime, these violations are not important. 

The scattering matrix 5* for this single channel quasi- 
one-dimensional (Q1D) system is 



S 



hi hi 
hi hi 



Bagwell H has obtained 



in = 1 



'2fei 



1 . r Lml , ,-rii ' 



where, Y r . 



and hi = - 



-1 
(.Vi 



K fl 



2m, 



(E n -E), k n 



f)]Sin[^( y i 

2m e 



(4) 

(5) 
(6) 



(E — E n ), E„ 



, -3-, 7 is the strength of the delta potential placed 
at a distance yi from the center of the quantum wire of 
width w. We find that a similar expression exists for the 
transition amplitude from the propagating mode (n=l) 
to the nth evanescent mode, given by 
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When the impurity potential is positive it can only sup- 
port scattering states. However when the impurity po- 
tential is negative, it can also support bound states, apart 
from the scattering states. For each n we get a sub-band 
of scattering states (E as a function of k n ). Similarly we 
get a bound state for each n, that are solutions to y| 



1 



2 for 







(8) 



For n=l we get a true bound state. The bound state 
for n=2 may or may not be a true bound state. If the 
impurity potential is such that the solution to Eq.JSJl lie 
in the energy range where n=l channel is propagating, 
then this bound state for n = 2, is degenerate with n=l 
scattering state and it becomes a quasi-bound state. 
We have also calculated to find 



[p(E) - Po (E)] total = 



p 



hv r , 



dx Cos(2k p x + rjp) 



V — 

" ^ hv r , 



E 



(9) 
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Here v p = hk p /m. J2 P denotes sum over all propagating 
modes and ^ e denotes sum over all evanescent modes. 
The 1st term on the R.H.S. is basically integrated change 
in the LDOS in the leads. Since the delta function po- 
tential is a point impurity, the integrated LDOS in the 
leads extends from — oo to oc. One can do the integration 
to find dxCos{2k p x + r) p ) = irCos(?] p )5(k p ). So it is 
zero unless the quasi-bound state coincides with k p = 0. 
In the case of extended impurities one can see that this 
term gives an unimportant small contribution that does 
not change with energy. Also the carrier concentration in 
the leads is normally large enough to screen away a small 
oscillatory LDOS completely. So the relevant quantity 
that appears in FSR is 



P(E)- P0 (E) = j: £ 



I ^pe I 



(10) 



This is actually the integrated local DOS around the im- 
purity site and decaying away from the impurity site all 
the way up to ±oo. In the regime of single propagat- 
ing channel p can take only one value, i.e. p = 1 and 
e = 2, 3, ...oo. Thus we can independently calculate both 
sides of Eq. Q starting from 1st principles, where we 
do not have to throw away dispersive behavior or energy 
dependence of self energy. 

We first present below a discussion and definition of 
the WKB regime for a Q1D system, because it is an in- 
teresting subject on its own. When the incident electron 
propagates in a potential where the wavefunction changes 
very slowly in space then very little reflected wave is gen- 
erated and that is taken to be the WKB regime jig . 
So a delta function potential in one-dimension (ID) has 
a WKB regime at higher energies, when the reflection 
probability is very small. In the inset of Fig.l, where 
we plot |rn| 2 versus incident energy we find that there 
are three regimes. One is to the left of point Pi where 
fn| 2 is large and also strongly energy dependent. The 

H 2 k 2 

other is between the points Pi and Qi where ^— 1 >> 7. 
These two regimes can be seen in ID scattering (e.g., a 
delta function potential in ID) and are the non-WKB 
and WKB regimes, respectively. The third regime is to 
the right of the point Qi, where again \ fn\ is very small 
and is hence a WKB regime, but the energy dependence 
of |fn| 2 is very large. Such a regime cannot be seen 
in ID and is a specialty of Q1D. So the energies that 
lie to the left of Pi is the non-WKB regime, where the 
electron feels the potential very strongly and is almost 
entirely reflected back. Energies to the right of the point 
Pi correspond to the WKB limit. Although, the system 
considered here is a Q1D system, corresponding to a scat- 
terer in Q1D, there is an energy dependent scatterer in 
ID [ltA 1 2(1 l21j . The bound states and scattering states 
of these two potentials are identical and this is an exact 
correspondence, valid in all regimes, quantum or semi- 
classical. And so when the reflection probability is small 
in Q1D, it is also small in the corresponding ID poten- 
tial. Then all the notions and results of WKB regime 



that we are familiar with in ID are also true in Q1D. 

In Fig.l we find a large deviation of ir[p(E) — po(E)] 
(dotted curve) from (solid curve) at energies in the 
non-WKB regime (left of Pi). This is similar to what 
is seen in ID, 2D or 3D. In the WKB regime, that is to 
the right of the point Pi, although | f\\ | 2 is very small, 
its energy dependence is not as negligible as that of a 
potential in ID (eg, a delta function potential in ID or 
a square well in ID). Energy dependence of |fn| 2 auto- 
matically implies energy dependence of t p (or S a ), i.e., 
dispersive behavior. So there is an appreciable difference 
between it\p(E) — po{E)] and 

Let us now analytically analyse the curves (solid and 
dashed) to the right of Qi. In this region K2 — ► 0. From 
Ea. (|10(l we find that only the 1st term in the series is 
relevant. That is 



7T [p{E) - po{E)] Q diverges as 



2tt Tn 1 

hvi T 2 2 K2 



^0 

(11) 

Thus the strong energy dependence in the energy beyond 
Qi is due to the rapid population of the second subband 
through evanescent modes, as it approaches its propagat- 
ing threshold. Hence, unlike the Fano resonance, this is 
not a quantum interference effect. This is like the Van 
Hove singularity at the band edge. Similarly one can find 



dSj 



dE 



dE 



arg{ti 



(12) 



diverges identically. Note that although arg(tn) can 
have a discontinuity, the derivative exists at all energies. 
Essentially the right derivative and left derivative is the 
same at the discontinuity. Hence we prove FSR is exact 
as K2 — > 0. This is understood when we note that when 



K 2 — ► 0, fn I goes to zero at the band edge 



Also 



it is known that when ryy = then t p maximizes |12| . 
and energy dependence of gp 1 being negligible, ^§r=0 . 
Thus all the deviating terms being zero, the FSR is valid 
around the diverging DOS at the band edge. 

For strong negative potentials, such that the bound 
state for n — 2 is below the propagating threshold of 
n = 1, the curves look similar to that in Fig. 1. For 
the negative 8 function potential with the bound state 
for n = 2, in the propagating regime of n = 1, we have 
plotted the two sides of Eq.QJ in Fig. 2. |iii| 2 is shown 
in the inset. Note that |tn| 2 shows that at the point P 
the system is in extreme non-WKB regime where \tn\ 2 
goes to zero. At this energy there is a quasi-bound state 
and there is strong energy dependence of scattering ma- 
trix elements as well as self energy. According to earlier 
stated results pj , there should be violation of FSR here. 
The peak in tt[p(E) — p (E)} at P occurs due to the quasi- 
bound state. We also see that at this very point P there 
is an exact agreement between L.H.S. & R.H.S. of Eq.JIJ. 
This can be even verified analytically. Substituting the 
bound state condition given in Eq.(|SJ) into -j£ as well as 
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in ir[p(E) — po{E)] separately, we get 

This agreement between and ir[p(E) — po(EY\was 
argued to be equal for the case of a stub in ref at 

the transmission zero where the _?Sr term in Eq. J2J) was 
dropped from the very beginning 6] . Dropping the energy 
dependence of S a in non-WKB regime and verifying the 
validity of FSR is rather meaningless, as the violations 
do come from the energy dependence of S a . Even af- 
ter including these terms we get exact agreement at the 
transmission zero for the negative 5 function potential 
in a quantum wire, although it is in extreme non-WKB 
regime. The reasons are as follows. At the transmission 
zero, since there is a quasi-bound state, E° becomes min- 
imum and ^g- = 0. All these arguments are also true 
for the stub. 

The multichannel case was analysed in Ref.^J. When 
two modes are propagating, the bound state coming from 
the 3rd sub-band, can be degenerate with the propagat- 
ing channels. At that point there is no agreement be- 
tween ^ and n[p(E)-p (E)}. While n[p(E) - p (E)} is 
not only positive definite, it also has a sharp peak mean- 
ing enhanced DOS at the energy corresponding to the 
bound state. The peak is completely missing from 

and turns out to be negative. Thus the disagreement 
is not just a quantitative one, but is a qualitative one. 
Thus the purpose of this work was to verify FSR in sin- 



gle channel quantum wires in the presence or absence of 
Fano resonances. Fano resonance is a very general feature 
of quantum wires. At the Fano resonance all the quanti- 
ties are strongly wave vector dependent as it is a purely 
quantum interference effect. Never the less, FSR is exact 
at the Fano resonance. This is contrary to the known fact 
that FSR is valid in semiclassical regimes where there is 
no strong dependence on wave vector. The exact agree- 
ment of the FSR in spite of the strong wave vector de- 
pendence is due to the fact that at the Fano resonance 
there is a quasi bound state that leads to a minimum in 
the self energy. Away from this quasi bound state there 
are strong violations. These are true for any negative 
potential in Q1D and the potential considered here and 
the associated calculations make this clear. For positive 
as well as negative delta function potentials, there is also 
strong wave vector dependence, close to the upper band 
edge of single channel propagation. This is due to the 
rapid population of the first evanescent mode at its prop- 
agation threshold and does not depend on the existence 
of Fano resonance. dY,/dE=0 here because of the perfect 
transmission at the band edge, and hence the agreement 
in FSR. So the former case of agreement in the peak is 
an agreement in purely quantum regime, while that in 
the case of the latter peak is in the semiclassical regime. 
Away from the peaks there is always violation. It may 
be interesting to work out some extended potentials in 
Q1D P. 
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Figure Captions 
Fig. 1. The dashed curve gives ir(p — po) and the solid 
curve gives -5i In Det[S]). Both the functions are 

plotted versus EW 2 using yi = .21W and 7=1. In 
the inset the corresponding |rn| 2 is plotted. We have 



considered 500 evanescent modes. 

Fig. 2. The solid curve gives w(p — po) and the dashed 
curve gives .bilnDet[S\). Both the functions arc 

plotted versus EW 2 using yi — .21W and 7 = —1.5. For 
this value of 7 there is a quasi-bound state at EW 2 = 
36.1022. In the inset the corresponding \tn\ 2 is plotted. 
We have considered 500 evanescent modes. 



